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We propose a new motivation for the stability of dark matter (DM). We suggest that tlie same non- 
abelian discrete flavor symmetry which accounts for the observed pattern of neutrino oscillations, 
spontaneously breaks to a Z2 subgroup which renders DM stable. The simplest scheme leads to 
a scalar doublet DM potentially detectable in nuclear recoil experiments, inverse neutrino mass 
hierarchy, hence a neutrinoless double beta decay rate accessible to upcoming searches, while ^13 = 
gives no CP violation in neutrino oscillations. 
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Introduction The existence of dark matter (DM) 
plays a central role in the modehng of structure formation 
and galaxy evolution, affecting also the cosmic microwave 
background. Despite the strong evidence in favor of DM, 
its detailed nature remains rather elusive. Viable particle 
physics candidates for dark matter must be electrically 
neutral, and provide the correct relic abundance. There- 
fore they must be stable over cosmological time scales. A 
simple way to justify the stability of the DM is by assum- 
ing some parity symmetry Z2, which might arise from the 
spontaneous breaking of an abelian U(l) gauge symme- 
try [l|-y| ^, or from a non-abehan discrete symmetry, as 
might be the case in some string models [J]. 

Non abelian discrete symmetries are motivated by 
neutrino oscillation data [5|, |6[ . Here we propose that the 
same symmetry explaining neutrino mixing angles is also 
responsible for the dark matter stability. In our simplest 
type-I seesaw [7| realization the flavor symmetry A4 
spontaneously breaks to Z2 providing a stable DM 
candidate. We extend the scalar sector of the standard 
model by adding three Higgs doublets transforming as a 
triplet of A4 we show that there is a consistent pattern 
of vacuum expectation values (vcvs) for which only one 
of the three extra Higgs doublets takes a vev, while 
the other two give rise to the dark matter candidate. 
The model accounts for the observed pattern of mixing 
angles [8| indicated by current neutrino oscillation data, 
predicting ^13 = and inverted spectrum of neutrino 
masses. It will therefore be tested in upcoming double 
beta and neutrino oscillation searches [9| , while the dark 
matter has potentially detectable rates within reach of 
nuclear recoil experiments. 

Model We assign matter fields to irreducible repre- 
sentations of A4 , the group of even permutations of four 
objects, isomorphic to the symmetry group of the tetra- 
hedron. All elements are generated from two elements S 



and T with S"^ = T^ = {STf = 1. A4 has four irre- 
ducible representations, three singlets 1, 1' and 1" and 
one triplet. In the basis where S is real diagonal. 
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one has the following triplet multiplication rules, 

{ab)i = ai^i + 0262 + 03^3 ; 

{ab)i' ~ ai6i +0^0262 +w^a363 ; 

{ab)i" = ai6i + ^^0262 + ^0363 ; 

(a6)3i = (02^3,0361,0162); 

(06)33 = (0362,0163,0261), 



(1) 



(2) 



where oj^ = 1, a ^ (ai, 02, 03) and b ~ (61, 62, 63). We as- 
sign the standard model Higgs doublet H, to the singlet 
1. and we assume three additional Higgs doublets trans- 
forming as an A4 triplet, namely 77 = (771, ?72, '73) ^ 3. 
We have four right-handed neutrinos, three transform- 
ing as an A4 triplet Nt = (iVi, A^2, -^3), and one singlet 
A'4. The lepton and Higgs assignments of our model is 
in table|Tl The resulting Yukawa Lagrangian is 
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^ In supersymmetry a viable DM particle is the neutralino, whose 
stability stems from the imposition of the so-called R-parity. 



TABLE I: Summary of relevant model quantum numbers 



C = yeL,llH + y^L^llH + yrLrl'^H + 

+y'(UXNTv)i + y!^L^{NTv)i" + ylLr{NTv)v + 
+ylLeN4H + MiNtNt + M2N4N4 + h.c. 

This way H is responsible for quark and charged lep- 
ton masses, the latter automatically diagonal. Note that 
we do not discuss the quark sector, assumed to be blind 
to A4, namely all left and right-handed up and down- 
type quarks transform trivially under A4 , their mass and 



mixing hierarchies might arise from an extra family sym- 
metry, for example, Frogatt- Nielsen-like [10|. Neutrino 
masses arise from H and 77, see below. The relevant terms 
of the scalar potential are of the form 



V 






^i]jH^ H + \i[H^ HY + X2[rihi]l- 



1" + \iWri^]v[m]i" + KW^iMwh 

+ A6(E.,,M^l3jw]3, + h.c) + XjWif^^H^H+ 

+ \'j[ri^H]H^ + \i{[i^^'n%HH + h.c) + 

+ A9 {Wvh^V^H + h.c) + A^ [Wih.V^H + h.c) + 

+ Aio {[i^hj%,i]H + h.c) + A'lo {Wv%^i^H + h.c) + 

+ All [ry^TytJi [7777] 1 

(3) 
where i,j = 1,2, and [...Js; means the product of two 
triplets contracted into a triplet of ^4, see eq. ^, [■■■]i 
means the product of two triplets contracted into a sin- 
glet of A4 and so on. Note that [77^?7]i,i'a" = [7777^]i,i',i", 
[mhi = [w]32 and so on. 

We have studied the minimization of the potential V 
solving the equations dV/dvi = where Vi are the vevs 
of the fields iJ, 771,772 and 773. For simplicity we assume 
real vevs. We have checked that for suitable parameter 
choices of the potential V , an allowed local minimum is 

(i7")=7;„^0, (77?) = 7;„^0 (%%) = 0, (4) 

with the eigenvalues of the Hessian d^V/dvidvj all posi- 
tive. 

Note that the alignment (77) ^ (1,0,0) breaks sponta- 
neously A4 to Z2 since (1,0,0) is invariant under the S 
generator in eq. ^. The Z2 is defined as 



(5) 



This residual symmetry is responsible for the stability 
of our DM candidate and the stability of the minimum. 
Note that the potential cannot break spontaneously A4 
into Z3 because in this case the alignment (7/) ~ (1, 1, 1) 
is not a minimum unless a fine tuning in the parameters 
A9+A10 = is assumed. This atractive feature reminds of 
the inert doublet model [ll| , with the difference that here 
it follows naturally from the underlying flavor symmetry 
which accounts for neutrino oscillations. 

We have four Higgs doublets^ giving three physical 
charged scalar bosons, plus four neutral scalars, and three 
pseudoscalars. After electroweak symmetry breaking we 
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^ Lepton flavor violating processes are suppressed by the large 
right-handed neutrino scale. 



can write 
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(6) 



There are 3 physical charged scalar bosons, 4 CP even 
and 3 CP odd neutral scalars. The mass of the neutral 
scalar fields is block diagonal with the standard model 
Higgs h mixed with hi, but not with the scalar fields 
with zero vev's /i2,3- 

Dark matter The lightest combination of the stable 
scalar fields ft.2 , ^3 plays the role of our dark matter par- 
ticle, which we will denote generically by rjohi. We list 
below all interactions of 77/) j\/ : 



1. Yukawa interactions 

Vdm ^1^2,3 , 

where i = e, fi, r. 

2. Higgs- Vector boson couplings 

VomVomZZ, VdmVomWW, 
VnM^2,3Z . 

3. Scalar interactions from the Higgs potential: 



(7) 



(8) 



Vdm AiA2hi , 

Vdm ^2^3/13 : 



VomV 



DAI I DM 



hh , 



Vdm AiA3hi , 
77^,, AiAaft. , 
Vdm hih3h 

VdmVdm l^if^i 



(9) 



After electroweak symmetry breaking, the vevs Vh and 
Vri are generated, so that additional terms are obtained 
from those in Eq. ([9]) by replacing h — > Vh and hi — > W;,. 
The flavor symmetry A4 is broken down to the residual 
Z2 symmetry in Eq. ([5|), implying the stability of our 
dark matter candidate. As we will see, despite the 
many mass and coupling parameters appearing in the 
potential, eq. ([3]), for A/^ ^ Mz, only two determine 
the relic dark matter abundance and its direct detection 
rates. 

Relic Density Assuming that our DM candidate 
arises as thermal relic in the early universe, one of 
the most important requirements one must check is its 
relic abundance. For definiteness we require that rjjjM 
makes up all the observed DM. For M^ ^ M^ the most 
important annihilation and coanihilation processes are 
those with vector bosons, though for large A > 172, where 
16A2 = (A7 + A'7 + 2Xs)^ + (2A2 - A3 - 2A4 + A^ -f 2A5 + 
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FIG. 1: Feynman diagrams relevant for direct DM detection. 
Elastic scattering (left) is generically more important than 
inelastic (right). 



A5 + Ae)^, annihilation into Higgs bosons plays an im- 



portant role, see En 
approximated as [12l | 

n-DMJT) 
s{T) 



(|TT|) . The DM abundance can be 



180 



n g^ Afpi Tf{(TAv) 



(10) 



where 



A/„ 



26 and g, = 106.75 + n is the number of SM 
degrees of freedom plus 1 < n < 12 degrees of freedom 
arising from the extra scalars, and Mpi = 1.22 x lO^^GeV 
is the Planck scale. The cross section for tjjjmVdm -^ 
VV where V are vector bosons in the limit of massless 
final states, is given by [12| 



(CAW) 



3g: 



'4 



Qgigi.Y'+Ay 



2567r M2 



(11) 



Y 



where 

92 = \/47ra/(l - A/^./A/|) and gy 



1/2 is the weak hypercharge, 
^A-KaMz/Mw 
From these equations it follows that, in order to provide 
the correct relic abundance Vl-oyih? = 0.110 ± 0.006 i.e. 
tidm/s = (0.40±0.02)ey/M^ [3], a correlation between 
the mass of the dark matter M^ and the quartic coupling 
constant A is required. For simplicity if we take the limit 
of small A we obtain a mass for the DM candidate of 
A/,, « 0.51 TeV. For large A values we have that the DM 
mass Mrj scales as A. 

Direct detection The quartic couplings rj^rjH'^ H and 
rj^rj'' HH give an interaction of the DM candidate with 
the nucleon through the interchange of the SM Higgs 
boson. Hence our DM candidate can be detected through 
the elastic scattering with a nucleus tjdmN — > r]oMN via 
the exchange of a Higgs, or through inelastic scattering 
with a nucleus tjdmN -^ AN with the exchange of a Z 
boson, see Fig. [I] where A is the lightest pseudoscalar, 
in general a mixture of A2 and ^3 . 

Barring fine-tuned choices of parameters for which the 
threshold for inelastic scattering opens up, the detection 
will be dominated by the elastic process, whose cross sec- 



tion is given by [14| 

(Tel (nucleon) w A 



1 + (tan/ 



lOOGcVV 
Mh J 



/ 50GeV y / 



5 X 10"*^ cm^ , 



(12) 



where tan/3 = Vh/vrj. Note that all uncertainties as- 
sociated with the nuclear form factor in Eq. (fT2|) . have 
been neglected. From the requirement of correctly re- 
producing the relic density, Eqs. (fTU]) and pT|) . one can 
find an expression for A as function of the DM mass, 
M,f. Using this relation and eq. (fT2|l one can plot the 
estimated cross section for the direct detection for each 
value of tan/3 and mass of the Higgs, Mh, as illustrated 
in Fig.[2]'^. The figure has been generated using [15| and 
compares the experimental sensitivities with our model 
expectations, fixing mn = 120 GeV and three tan/3 val- 
ues. This choice of Higgs mass is motivated by the LEP 
bounds rriH > 114GcV, which however is not strictly 
valid in our model due to the additional Higgs doublets. 



DATA listed top to bottom on plot 
XENON102007,iii««re(iLet[ from Xe cube 
CDMS: Soudan 2004-2009 Ce 
ZEPLlNI(vt 3,widi PMT iipaide) Pioj. Sens. 

'"■■ ■ d sensitivity: fflOOkt-d 5-30 keV.45?r elf. 

SiiperCDMS'-KlOkgatSNOLAB ' 
XENONIOO upgrade mjected sensitivity: 60,000 kt-d 5-30 keV, 454 efl. 
^ IT projectedsensidvity: 3 ton-yr, 2-30 keV, 454 eff. 
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FIG. 2: Elastic DM scattering cross section with a nucleon 
versus DM mass. We compare present [Ig, [13] and future 
[la . Il9l | sensitivities with our model expectations, for niH = 
120 GeV and tan/3 = 0.5, 1, 5 (grey sohd lines). 



Neutrino phenomenology Our model has four heavy 
right-handed neutrinos, and is a special case, called (3,4), 



^ Here we focus on the region Mr] ^ Mz 
light DM will be treated elsewhere. 



The interesting case of 



of the general type-I seesaw mechanism (20| . After elec- 
troweak symmetry breaking, it is characterized by Dirac 
and Major ana mass terms given as 

xi j/i 
ruD^ \ X2 I ,MR = diag{Mi,Mi,Mi,M2) 
X3 

(13) 
so that the hght neutrinos get Majorana mass by means 
of the type-I seesaw relation rrii, = —rriDs^^M^^ 4'^i)3x4 
the light-neutrinos mass matrix M,^ being given as 
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It falls within the class of scaling matrices introduced in 



Ref . j21| . This form of the light neutrino mass matrix 
has an inverse hierarchical neutrino mass spectrum and 
a zero eigenvalue with TO3 ~ and corresponding eigen- 
vector (0, —e/h, 1)^ implying a vanishing reactor mix- 
ing angle 6*13 = 0. One can see explicitly that the solar 
and atmospheric square mass differences and mixing an- 
gles indicated by neutrino oscillation data [8| can indeed 
be fitted by taking, as an example, the tri-bimaximal 
(TBM) ansatz [22|. When b^ e and y^ ^ 2c^ - ae the 
neutrino mass matrix Eq. (|14p is /z — r invariant yield- 
ing maximal atmospheric mixing, sin ^23 = 1/2 and 
M^4i + M^,i3 = A/^,22 + A^^,23, which gives the TBM 
value of the solar angle, sin 0^2 = 1/3, in good agreement 
with experimental data within one a. The eigenvalues 
are {m-i, ?Ti2, w.3} = {2ac + 2c^, 2c^ — ae, 0} , which can 
fit the two mass-squared differences required to account 
for the observed pattern of neutrino oscillations. By re- 
laxing the condition 6 = c and y^ = 2c^ — ac one generates 
deviations from the TBM limit, while keeping 6*13 = 0. 
Note the above imples a neutrinoless double beta decay 
effective mass parameter in the range 0.03 to 0.05 cV at 
3 a, within reach of upcoming experiments |23| . 

Conclusions In summary we have suggested that 
DM stability follows from the same non-abelian discrete 
fiavor symmetry which accounts for the observed pattern 
of neutrino oscillations. In the realization we have given 
we have an A4 symmetry which spontaneously breaks to 
a Z2 parity that stabilizes a scalar doublet dark matter, 
potentially detectable in nuclear recoil experiments, 
as well as accelerators. Despite the complexity of the 
scalar potential, in the heavy dark matter limit both 
the relic dark matter abundance and its direct detection 
cross section depend just on the DM mass and a single 
coupling strength parameter. The model is also mani- 
festly unifiable and agrees with electroweak searches as 
well as precision tests, as will be shown elsewhere. Our 
simple example gives Oi^/3/3 rates accessible to upcoming 



experiments and no CP violation in neutrino oscillations. 
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